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Abstract: In this paper, data compression performance of high SLFOR fixed windows based cosine 
modulated filter banks is presented. Three fixed windows with high SLFOR (decay^reapr^an or equal to 
-24 dB/octave) are selected for this study. The designed prototype filters are optimi: 
banks are developed. These filter banks are used for compression of ECG signal data 
is a dictionary based scheme, is used for sub band coding of the signal. 
computed for quantitative assessment of reconstructed ECG signal. \J 
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I. Introduction 




ismn Ratio, LZW algorithm. 



Electrocardiogram (ECG) signal represents the electrical \jjT\j heart and is vital in several clinical 
diagnoses. The need for ECG data compression arisrt|ui various applications like ambulatory ECG 
monitoring, patient databases in hospitals, medical educa&Jr systems, ECG telemetry, etc. [1]. The digitized 
ECG generates massive volume of data, henodLrailf^Lre* an effective compression technique scheme for 
transmission/ storage. The main goal of anv Amp^ession technique is to achieve maximum data reduction 



while preserving the significant sign 
approaches for ECG data compression: 
methods. Direct methods like TP, AZTE 1 
which utilizes correlation of neighb 
compression methods apply a transfi 
strategy is to minimize the 
coefficients, compactly repr 





features upon reconstruction. There are three 
ethods, transform methods and parameter extraction 
OfcRTES, FAN, SAPA, etc. are sometime domain based schemes 
mples to reduce the redundancy [2-4]. The transformational 
o analyze the energy distribution of the signal. The compression 
%f coefficients, which can represent the input signal. These selected 
the input ECG signal, are transmitted/ stored, while other redundant 



samples are discarded. So*fee oWthe transformational compression methods are Fourier transform, discrete 
cosine transform (DCT)^Alsh transform, Karhunen- Loeve transform (KLT), wavelet transform and filter 
banks [5]. The paramerfc^extraction methods are prediction based coding, where near prediction behavior 
of ECG signals areiusVd to develop subset auto regression (SAR) model for long term prediction [6]. 



1 /-\ rjn)/^ . 



FJn] 














m 


w 




U4 





Fig. 1. M-Channel maximally decimated filter bank. 
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Multirate filter banks (Fig. 1) find variety of applications in sub band coding, trans multiplexing, image and 
video or audio compression, spectral estimation, bio signal processing, and adaptive signal processing. The 
fundamental block in implementation of such applications is cosine modulated filter banks (CMFBs), where 
analysis and synthesis filter banks are derived by cosine modulating the low pass prototype filter. Hence the 
entire design of the filter banks is reduced to design the prototype filter. Therefore, during the design 
phase, it is required to optimize the coefficients of prototype filter only. This significantly reduces the 
complexities and computational overheads. 

To obtain high quality reconstruction, the linear phase low pass prototype filter P(e0) must satisfy 
following two conditions as much as possible [7] -[8]: 



p(0 



= 0J&r|fli| > irjM 
= 1 for Q<(o<7rM: 



(1) 



where T(e J ) = I 
Jt=0 



2M-1 




M is the number of channels. 

If (1) is satisfied exactly, there is no aliasing between ^jpnadjacent bands, while if (2) is satisfied, all 
amplitude distortion is eliminated in the combi^flWkalfsis /synthesis filter bank system. Aliasing between 
the adjacent bands is eliminated by selecting^ppiTOrJate phase factor in the modulation. Hence the filter 
bank that provide approximate or neaj^peWc^econstruction (NPR), is designed, that approximately 
satisfies the constraints laid down in (1) ad (1^ Design approaches involve linear optimization [7] -[9] as 
well as nonlinear optimization [10] of diff^reiwobjective functions. 



The prototype filter p(n) arj 
side lobe fall-off rate (SLp6r" 




Filter Bank Design 

ed with windowing technique, using high SLFOR fixed windows. High 
a spectral parameter, which is important in many applications. In beam 



forming applications, a SLFOR ensures a better rejection of far end interference [11]. Superior far end 
Dafc% eri 



attenuation for sto 
[12]. In speech 

[13]- 



r stop* b, 



energy may be achieved in filters, designed with high SLFOR window function 
nal processing applications, it minimizes the energy leaks from one band to another 



Nutgall [14] proposed some 3-and 4 non-zero term windows with high SLFOR characteristics. These 
windows exhibit rapid decay in response and have some higher derivative continuous for all t. These 
windows w(n) are defined for the DFT by: 



w{n) = ci -flj cos( — nj + a^ cos( — 2n)-fl 3 cos( — 3^); 

N N N 



with n = 7 l 7 ...,N-l 



(3) 



The coefficients of these windows are given in Table-I. 
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Table I. Window Parameters and Coefficients 



Window 


Max 
Sidelobe 
attenuation 


SLFOR 

dB/ 
octave 


Coefficients 


4-temi with, continuous 


-60.95 


-42 


ag= 10/32, 3.1 


fifth derivative 






=15/32., a 2 = 


(4TC5D) 






6/31 a 3 =l/32 


4-temi with continuous 


-82.00 


-30 


ao= 0.338946. 


third derivative 






a 3 = 0.481973, 


(4TC3D) 






a 2 =0.1G1054, 
3^=0.018027 


3 -term with continuous 


-4G.74 


-30 


a = 0.375. 3,= 


third derivative 






0.5, 3] = 0.125, 


(3TC3D) 






3 3 =0.00 



For prototype filter p(n) =hid(n)-w(n) , where hid(n) is the shifted impulse response, 
Cutoff frequency w c is (w s + w p )/2 and transition bandwidth is (w s - w p ) are calculated. 



Lilatec 



d.> 



low pass filter. 



The cutoff frequency w c of low pass prototype filter is taken as the variable parameter, which is optimized in 
such a way so as to yield the minimum value of objective function. The !|£rrfl^Kjective function as in [7] is 
used and the algorithm evaluates w c such that the 3-dB cutoff poir£ of opwpized prototype filter is located 
approximately at 2M. j£ • ' 

2 



F - max 



The golden section search method [15] is 
frequency w c of prototype filter [16] . 




iJimize the objective function by adjusting the cutoff 



The impulse responses of the analysis mhemBanks hk(n) and synthesis filter banks fk(n) are the cosine 
modulated versions of the optimized pmtaffpe filter p(n) [17H19]. 
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Q<n<N-l, Q<k<M-l (5) 
Fig 2. Frequency response of optimized 32-channel cosine modulated 
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analysis filter bank, with prototype filter designed using 3-term with continuous third derivative (3TC3D) 
window. Only first eight bands are shown to highlight the high fall-off rate of the designed filter banks. 

III. Sub Band Coding Scheme 

A subband coder splits the input signal into a collection of approximately disjoint frequency bands. If the 
resulting subbands have the same extent in the frequency domain, the subband decomposition is said to be 
uniform. Since the bandwidth of each subband is reduced by an amount corresponding to the number of 
subbands, - say M, each subband, can be subsamples by a factor of M. Thus, the number of signal samples 
in the critically sampled subbands is the same as in the input signal. Since the critical information of 
various subbands is unequal, compression is obtained by representing (quantizing) the lesser informatic 
subband with a small number of bits. In case of ECG signals, normally, the signal energy is concentrated in 
the lower frequency subbands, implying that the higher frequency subband can be represented 
with a small number of bits or may even be discarded (Fig. 3). 



LZW (Lempel-Ziv- Welch) algorithm is a general compression algorithm capable of^^kmg on almost any 
type of data. It is generally fast in both compressing and decompressing data a^i^dts not require the use 
of floatingpoint operations. LZW is often referred to as a substitutional^ dieiionary-based encoding 
algorithm. The algorithm builds a 
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. 3. Subband codec system. 

data dictionary (also lmown^^ranslation table or string table) of data occurring in an uncompressed data 
stream. Patterns of data (substrings) are identified in the data stream and are matched to entries in the data 
dictionary. If the suljgtryStis not present in the dictionary, a code phrase is created based on the data 
content of the subs^rin^and it is stored in the dictionary. The phrase is then written to the compressed 
output stream _* A Y"IR a reoccurrence of a substring is identified in the data, the phrase of the substring 
already stored 1?N^e dictionary is written to the output. Because the phrase value has a physical size that is 
smaller than the substring it represents, data compression is achieved. 

Decoding LZW data is the reverse of encoding. The decompressor reads a code from the encoded data 
stream and adds the code to the data dictionary if it is not already there. The code is then translated into 
the string it represents and is written to the uncompressed output stream. The dictionary based 
compression schemes are normally faster than entropybased schemes, adaptive, and does not use any 
statistical model [20] -[21]. 



IV. Error Measures 

The quality of retrieved signal is measured using the Percentage Root Mean-square Difference (PRD), which 
is defined as [22] 
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PRD= l 1 ^^ xlOO (6) 

Where x[n] is the original ECG signal and y[n] is the reconstructed ECG signal and L its length. 

As PRD is heavily dependent on the mean value of signal, it is more appropriate to use the modified criteria: 



PRD1= \— xlOO 



I 
I 
pf=l 

Where x[n] is the mean value of the signal. 

If the value of PRDi is between o% to 9%, the quality c*P^(^o*istructed signal is either very good or good, 
where as if the value is greater than 9%, its quality group ta/mot be determined. In this paper, it is ensured 
that the PRDi value of reconstructed ECG signa^^l!Iltol to or less than 9%. 

The Signal to Noise Ratio (SNR) and PRDia^^^as [23]: 



SNR = - 2ologio (0.01PRD1) (8) 



Mean Square Error (MSE) is also on important parameter to evaluate the quality of reconstructed 

signal [24]. 
It is expressed as: 

S£ =-V| x(n) - y(n)\ 2 (9) 



is aiso onoK 



of the m 



The measure of the maximum differences gives an idea about the local effects. This can be achieved by 
means of the maximum amplitude error (MAX) and is defined as [22] 

MAX = max{ I x[n] - y[n]\] ( 1 0) 

=■; 

Compression ratio (CR) can be expressed as [25] 

The number of bits hi original 1-D ECG signal 

CR = (11) 

The total number of bits in the compressed data 

However, the clinical acceptability of the reconstructed signal should always be determined through visual 
inspections. 
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V. Design Examples and Results 

In first example, 32 channel cosine modulated filter banks are developed using optimized prototype filter. 
The prototype filter is designed with high SLFOR fixed window using windowing technique. The stopband 
frequency and passband frequency are chosen. Here is the roll off factor, influencing the overlap of 
subbands. Keeping the roll-off factor at 1.88, results in N= 380. Two ECG signal (records no. 100 and record 
no. 117) from ECG Arrythmia Database of MIT-BIH ECG database [25] are used for compression and 
decompression with LZW algorithm . These signals have been sampled at 360 Hz with 11 bits/sample. The 
1024 baseline is removed prior to subband coding of the signal. 



In second design, 16-channel cosine modulated filter banks are developed. The same Roll-off factor is at 
kept at 1.88, to obtain N= 190. In each case, a block of 1024 samples are used in the experiment. The results 
of both designs are given in Table - II. The original and reconstructed waveforms from filter banks designed 
with different high SLFOR fixed windows, are shown in figure 4. J^^^^ 

Results from Table-II show that Nuttalls rapidly decaying window (4TC5D window) orjgred superior signal 
compression performance than others. This is due to its high SLFOR of -^fe^iB Octave. However, the 
window has slightly lower sidelobe attenuation level of -60.95 dB, compared/fi^ 8^00 dB of 4TC3D. 



VI. Conclusion 



Performance of some high SLFOR fixed window functions bas^^iflyf* bank is evaluated for ECG signal 
compression. It is evident, that the sidelobe fall-off rate andsifceld&e* attenuation level are two important 
parameters that greatly influence the signal compression. AsQjifcussed, window with high SLFOR 

Table II. Signal Reconstruclfcj/r Error Measures 
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Original ECG, Reconstructed signal using: 4TC5D 
samples of ECG record no. 117 are shown) 



and high sidelobe attenuation level/^nVbr an excellent candidate in design of filter bank based data 
compression scheme. The fixed wilccp/s offer good performance in general but, has a drawback of 
inflexibility to provide desired sroeb^fid attenuation level. It can be concluded that high SLFOR property 
and good sidelobe attenuatirai are desirable in achieving a good signal compression ratio. 
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